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^ ■ Abstract 

Oh: 

•^ . We show that strong non~Hnear electron transport in the ferromagnetic insulating (FMI) state 

O^ I of manganites, responsible for phenomena such as colossal electroresistance and current induced 

resistance switching, can occur due to a hot electron effect. In the FMI state, which we show is 

I ■ an insulator with a Coulomb gap, the temperature of the electron and lattice baths can decouple 

U ■ 

c/5 , at high input power levels, leading to heating of the electron bath. Parameters of the hot electron 

+-* ■ 
J2 [ effect model were independently determined via time dependence experiments and are in good 

agreement with the experimental values. 
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Manganites Li„xAxMn03, where L and A are ions of a trivalent rare-earth and a divalent 
metal atom respectively, are known to exhibit a colossal magnetoresistance. Manganites 
also show electroresistance (ER) and related resistance state switching effects^. There is 
a resurgence of interest in these materials due to the possibility of their usage in memory 
devices^ based on electric (as opposed to magnetic) control of their electrical resistance. 
These effects are manifestations of a strong non-linear conduction (NLC)^ in these materials, 
wherein their resistance is a strong function of the bias current. NLC is generally observed 
in manganites having hole concentration 'x' such that the ground state is insulating, e.g., 
in the charge and orbitally ordered insulating ground state observed at x ~ 0.5, or in 
the ferromagnetic insulating (FMI) state, typically observed for x ~ 0.3. Though NLC 
phenomena for different ranges of 'x' are ostensibly similar, their physical origins are likely 
different. In this paper, we investigate NLC in the FMI state of manganites and suggest 
a mechanism that may underly it. NLC in the FMI state is responsible for a "colossal" 
electroresistance (CER) and related electric current induced resistance switching (CIRS}^»^. 
We show that NLC and electronic properties in the FMI state (an insulator with a Coulomb- 
gap) are intimately related. The NLC is modeled using hot electron effects, wherein exist 
separate temperature scales for the electron and the phonon/lattice baths. To our knowledge, 
this is the first proposal that such hot electron effects occur in manganites. 

We study single crystals of two widely different manganites, Lao.82Cao.i8Mn03 (LCM018) 
and Ndo.TPbo.sMnOs (NPMO30). As temperature T is decreased, the samples undergo 
paramagnetic-ferromagnetic transitions at T = Tc — 165 K and ^ 150 K respectively, and 
enter their FMI states at T = Tfmi — 100 K and ^ 130 K respectively (see left inset of 
figure [U the data were obtained at very low power bias). Examples of NLC data are shown 
in figure [TJ The strong dependence of resistivity p on current density j is evident. Extensive 
NLC data in these samples is available^>^. 

The transport in the FMI state is unique: it is activated below Tfmi; however, the 
FMI state is not a conventional band or Mott-Hubbard insulating state that have a gap 
in the density of states (DOS) near the Fermi level Ep. Low temperature heat capacity 
measurements of both LCM018 and NPMO30 reveal the presence of a finite electronic 
contribution 7T to the heat capacity. Such a situation can arise in a Shklovskii-Efros- 
insulator with a soft-gap (Coulomb-gap) Acq? and indeed, we find that p in the FMI state 
has a Shklovskii Efros variable range hopping (SE-VRH) T dependence: 
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p = p'^exp\-^\ (1) 

where, pg = poT^^^. Here, Tq represents an energy expressed in Kelvin. In right inset of 
figure [1] we show a plot of T~^/^ Inp versus T^^/^ for T ~ TpMi to show that the p follow 
equation [T] below TpMi- The experimentally estimated values of po. To, and Acg are listed 
in table [H The presented values of Acg are in agreement with direct estimation of Ace by 
angle resolved photoemission spectroscopy (ARPES)^. 

We propose that NLC effects in manganites arise due to heating of the electrons because 
of the presence of a thermal conductance, that is finite, between the electron bath and the 
phonon bath. In this situation the electron bath temperature Te will get decoupled from that 
of the phonon bath Tph when the power P input is large. Direct measurements of Tph have 
established that during measurements using higher P, Tph does not undergo any substantial 
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rise (~ 5 K) , i.e., the observed NLC cannot be due to heating of the sample because of Joule 
dissipation. Within the present model, shown schematically in figure Ej the NLC effects arise 
because of a heating of the electron bath. The energy exchange between the two baths is 
limited by an effective thermal conductance Ae_ph, modeled empirically as Ae_ph = AqT". 
Thermal conductances generally increase as the temperature is raised, implying a > 1. The 
phonon bath, in turn, is linked to the temperature controlled base through a finite thermal 
conductance. However, since it is known"^ that Tph remains close to the base temperature, 
we ignore the finite thermal resistance between the phonons and the base. The heat fiow is 
thus limited by Ac-ph- Using the relation, AT = (Tc — Tph) = P/Ac-ph, we obtain a working 
relation between Tph, Te, and P: 

Tph"^^ + ^P] (2) 

The values of Tq and po being known, at fixed Tph, equations [2] and [U together express 
p as a function of P in terms of only two parameters, namely a and Aq. Representative 
fits to this hot electron model for LCM018 and NPMO30 at indicated Tph are shown in 
figure [21 It can be seen that the simple model based on two parameters provide good fit 
to the data although there are deviations in the regimes of high power. The values of a 
and Aq based on the hot electron model are listed in table [B We discuss below that these 
parameter values can be related to experimentally determined numbers and are completely 



justifiable. In figure [2] we also plot on the right ordinate the variation of Tg as a function of 
P, as obtained via equation [21 It can be seen that at low power P < 10^^ W, Te is close to 
Tph as expected, but at higher power, Te increases rapidly and can be substantially greater 
than Tph for P > lO'^ W. 

The model proposed here is in the line of similar effects in heavily doped semiconductors^ 
with carrier concentration near the critical regime of metal-insulator transition. This regime 
is known to be dominated by localized electronic states and the conduction process is of the 
Shklovskii Efros type with a Acq in the DOS. The key difference between the two systems 
is in the scale of carrier concentrations and the scale of Acg- In doped semiconductors such 
as Si, these effects are manifested at carrier concentrations in the range of 10^^ to 10^*^ cm~^, 
and with Acg — 1 rneV. In the present hole doped manganites, the carrier concentrations 
are a few 10^^ cm~^, and Acq is about 2 orders of magnitude larger. Due to the high value of 
Acg these effects are visible in manganites at such high temperatures ~ 100 K, in contrast 
to heavily doped semiconductors where they typically occur below 1 K. 

The proposed model can be validated by an independent estimation of the two parameters 
a and Aq to compare with the values obtained by fitting the model to the NLC data. Reliable 
estimates of these parameters can be directly obtained in the following way: a relaxation 
time determines the energy transfer between the two baths (electron and phonon). In the 
temperature range we are working, the phonon bath heat capacity is much larger than that 
of the electronic bath. (For example, even at 50 K, the electronic heat capacity Ce is not 
more than 2% of the lattice heat capacity.) Thus the phonon bath acts like an infinite heat 
capacity bath to the electrons and in this situation, the relaxation time re_ph is given by: 

T-c-ph = -r-^ (3) 

^ic-ph 

Using Cc = 7T, one gets re_ph ~ ^T^^^ j Kq. 7 is known from the heat capacity mea- 
surements. To obtain a and Aq we need to know Xg.ph, which is measured by a step change 
experiment as described below. At a given temperature below TpMi, we increase bias cur- 
rent I by a step, causing a decoupling of Tph and Te. Since Te needs a finite time to reach 
equilibrium, the voltage V across the sample (which measures the resistance, which in turn 
depends on Te) lags behind the current step. We measure the time for V to relax to the new 
value. An example of such a voltage relaxation experiment for LCM018 is shown in figure [3l 



The voltage evolution follows a simple exponential function with a single relaxation time, 
which we argue is essentially the Te_ph: the fact that the Te changes with a finite relaxation 
time when a step power is applied to the electron bath would mean that the resistance, 
which is function of Te will evolve with time giving a finite response time for the voltage. 
(Note: there may be two other causes for the voltage to respond with finite time: (1) a 
capacitive effect, and (2) a finite thermal relaxation time for the lattice temperature to relax 
to the temperature controlled bath temperature, i.e., for the sample to thermalize. The first 
effect has been tested for and was found absent, and as for the second effect, reasonable 
estimates as well as measurements show that this relaxation is an order of magnitude faster 
than that seen in the present experiment.) Thus, the measured relaxation time is essentially 
Tc-ph, which was found to be have a weak temperature dependence. Since re_ph ~ t^"", the 
weak temperature dependence implies {a — 1) to be very small, i.e., a ~ 1. The numerical 
value of Te-ph — 200 ms. These facts together yield an estimate for Aq ~ 2 — 6 x 10"'^ J/K^s. 
The values of a and Aq determined from the fit to the data using the model (see table [T| 
are in good agreement with these estimates. Similar agreement was obtained for NPMO30 
as well. Thus the simple model is validated by an independent experimental estimation of 
parameters. 

We note that the electron-phonon energy relaxation which is parameterized through a 
single Te_ph, needs a finite Ce, i.e., a non-zero 7. This can only happen when the insulating 
state is created by localization type phenomena, i.e., when the sample is not a band insulator. 
Presumably therefore, the type of hot electron phenomena that we are proposing would need 
a specific type of insulating state, and may not be applicable in other types of insulators, at 
least not in the present form. 

The deviation between the simple hot electron model and the experimental data at some- 
what high power levels can be due to the presence of electronic inhomogeneities in the FMI 
state. Tacit in the simple model analysis is an assumption of uniform dissipation of power 
within the sample. This assumption is not strictly true for manganites, especially those 
having compositions close to M-I phase transition boundaries, as is the case for the com- 
positions studied here. The inhomogeneous conductivity can lead to inhomogeneous power 
dissipation within the electron system, and thus to an inhomogeneous distribution of Te. 

To summarize, we have investigated non-linear conduction in hole doped manganites 
in the hole concentration region where they show a FMI state. The insulating state in 



the presence of ferromagnetic order in the hole doped manganites Lao.82Cao.i8Mn03 and 
Ndo.rPbo.sMnOs is caused by electron localization due to Coulomb interaction and the charge 
transport is of the SE-VRH type. The value of the Coulomb-gap obtained from the data 
are large ~ 10^ meV. We find that a simple model of electron heating can explain the data. 
The increase in bias current densities, i.e., increase in input power, leads to a heating up 
of the electron bath which is coupled only weakly to the phonon bath via a finite thermal 
conductance. The parameters of the model, estimated through independent experiments, 
were found to be close to those obtained by fitting the model to the data. 
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Composition pQ Tq Acg ^o « 

(SlcmK-i/2) (K) (meV) (W/K^+^m) 
LCM018 5.6 X 10^6 2.5 x 10^ 200 2 x lO"'^ 1 
NPMO30 1.8 X 10-2 3.9 x 10^ 40 2 x 10"^ 1.3 

TABLE I: Shklovskii-Efros transport, pQ, Tq, and Coulomb gap AcG) ^-^d hot electron model, Aq 
and a parameters. 
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List of Figure Captions 

FIG.l: Representative resistivity p vs current density j data at indicated fixed phonon tem- 
perature Tph. Left inset: Temperature T dependent p. Tlie respective Tfmi's are 
indicated by arrows. Riglit inset: Sliklovskii Efros variable range flopping transport 
in tfie FMf state. 

FIG. 2: Experiment aliy observed variation of resistivity p as a function of power P for (a) 
LCM018, and (b) NPMO30, at indicated fixed plionon temperature Tph. Tlie re- 
spective fits to tlie hot electron model (shown inset) are also shown. Right ordinate: 
calculated variation (using equation [2]) of electron temperature Tg as a function of P. 

FIG. 3: A typical voltage V relaxation profile for LCM018 upon the application of a positive 
step change in current bias I at time = 10 s, at indicated phonon temperature Tph. 
The solid line is a fitted exponential decay curve used to extract the electron-phonon 
relaxation time constant Tc_ph. Inset: Tph dependence of re_ph. 
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FIG. 1: Representative resistivity p vs current density j data at indicated fixed phonon temperature 
Tph- Left inset: Temperature T dependent p. The respective TpMi's are indicated by arrows. Right 
inset: Shklovskii Efros variable range hopping transport in the FMI state. 
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FIG. 2: Experimentally observed variation of resistivity p as a function of power P for (a) LCM018, 
and (b) NPMO30, at indicated fixed phonon temperature Tph- The respective fits to the hot 
electron model (shown inset) are also shown. Right ordinate: calculated variation (using equation[2]) 
of electron temperature Tg as a function of P. 
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FIG. 3: A typical voltage V relaxation profile for LCM018 upon the application of a positive step 
change in current bias I at time = 10 s, at indicated phonon temperature Tph- The solid line is a 
fitted exponential decay curve used to extract the electron-phonon relaxation time constant rc_ph- 
Inset: Tph dependence of re_ph- 
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